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ABSTRACT 

The abundance of the most massive objects in the Universe at different epochs is a 
very sensitive probe of the cosmic background evolution and of the growth history of density 
perturbations, and could provide a powerful tool to distinguish between a cosmological 
constant and a dynamical dark energy field. In particular, the recent detection of very massive 
clusters of galaxies at high redshifts has attracted significant interest as a possible indication 
of a failure of the standard ACDM model. Several attempts have been made in order to 
explain such detections in the context of non-Gaussian scenarios or interacting dark energy 
models, showing that both these alternative cosmologies predict an enhanced number density 
of massive clusters at high redshifts, possibly alleviating the tension. However, all the models 
proposed so far also overpredict the abundance of massive clusters at the present epoch, and 
are therefore in contrast with observational bounds on the low-redshift halo mass function. 
In this paper we present for the first time a new class of interacting dark energy models 
that simultaneously account for an enhanced number density of massive clusters at high 
redshifts and for both the standard cluster abundance at the present time and the standard 
power spectrum normalization at CMB. The key feature of this new class of models is the 
"bounce" of the dark energy scalar field on the cosmological constant barrier at relatively 
recent epochs. We present the background and linear perturbations evolution of the model, 
showing that the standard amplitude of density perturbations is recovered both at CMB and 
at the present time, and we demonstrate by means of large N-body simulations that our 
scenario predicts an enhanced number of massive clusters at high redshifts without affecting 
the present halo abundance. Such behavior could not arise in non-Gaussian models, and is 
therefore a characteristic feature of the bouncing coupled dark energy scenario. 
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1 INTRODUCTION 

The origin of the observed acceleration of t he cosmic expan- 
sion JRiessetijJll998l : |Perlmutteret"Zlll999l : lAstier et aLlbood : 
iKowalski et al.l 2008) represents one of the most relevant open 
questions in astrophysics and cosmology. Besides the standard 
ACDM cosmological model, that assumes a perfectly homoge- 
neous and static energy density represented by a cosmological 
constant A as the source of the acceleration, alternative scenarios 
based on a dynamical dark energy (DE) field such as Quintessence 
( fWetterich 1988; Ratra & Peebles 1988; Ferreira& Joyce 199^ 
and k-essence jArmendariz-Picon et alj200Ctl200l[ ) , or on a modi- 
fication of General Relativity a t large scales (see e.g. IHu & Sawickil 
uOOTIuApplebv & Wellerl2010r) have been proposed in recent years. 

In order to efficiently distinguish among these different possi- 
bilities it is therefore crucial to devise observational tests capable 
of highlighting a possible dynamical nature of the DE field, thereby 
falsifying the cosmological constant picture. While present obser- 



vational constraints on the background evolution of the Universe 
seem to be fully consistent with a static behavior of the DE density, 
as expected for a standard cosmological constant, complementary 
and more sensitive tests of a possible time evolution of the DE field 
might arise from the study of structure formation processes, and in 
particular from the time evolution of density perturbations at dif- 
ferent scales. 



In this conte xt, the detection of unexpectedly massive clusters 

a'l]|2006l: IJeeet all 



at high red shift JMuUis et al.ll2005l; 



Brodwin et al. 



Bremer et i 
20101 : 



IJeeet all 1201 II : 



2009: Ros ati et al. 
Foley etaljT201lh 



has recently attracted significant interest as a 
possible indication of a failure of the standard ACDM cosmolog- 
ical model. In fact, starting from a Gaussian field of primordial 
scalar perturbations with an amplitude consistent with the most re- 
cent constraints from obs ervations of the Cosm ic Microwave Back- 
ground (CMB) radiation fKomat su et al.ll201 ih . the growth of den- 
sity fluctuations for a ACDM cosmology predicts an evolution 
of the halo number density that would imply an extremely low 
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probability for the detection of sucli clusters. In particular, for the 
mo st extreme case represented by the XMMU J2235. 3-2557 clus- 
ter jjee et al.ll2009l : iRosati et al.ll2009l) with an estimated mass of 
M324 = (6.4 ± 1.2) X 1O"M0 at ~ 1.4, a 2 - 3cr discrep- 
ancy with respect to the fi ducial ACDM model has been claimed 
by several authors (see e.g. Jee et al. 2009; Jimenez & Verde 2009; 
[HoIz & Perlmutter 2010; Hovle et al. 2011; Jee et al. 2011), al- 
though more conservativ e interpretations have b een provided by 
other analysis (as e.g. bv iMortonson et alj|201 ih . The main diffi- 
culties in driving sufficiently robust conclusions on a possible ten- 
sion with the standard ACDM model from such detections arise as 
a consequence of the observational uncertainties on the mass deter- 
mination of the clusters and on the estimation of the statistical sig- 
nifica nce of the cosmologica l volumes covered by t he cluster sur- 
veys JSheth & Diaferioll201 ll ; IWaizmann eTai]|201 ih . Such signifi- 
cance will need to be better clarified before claiming a discrepancy 
with respect to the predictions of the standard model. Nevertheless, 
it is in any case interesting to investigate whether such possible dis- 
crepancy could be alleviated by alternative cosmological scenarios. 
Several attempts have been made in this direction in the 
recent past. On one side, a deviation from statistical Gaus- 
sianity in the primordial density field has been shown to give 
rise to a higher number density of massi ve halos as compared 
to the standard Gaussian case (see e.g. Matarrese et al. 200(J; 
Grossietal. 2007; Jimenez & Verde 2009; Holz & Perlmutter 
2010; Cayon et al. 2010; ISartoris et ahlboiol ; iHovle et alj|201ll ; 
Lo Verde & Smith .201 1.1 . However, in order to significantly in- 



crease the probability of detection this approach requires a level 
of primordial non-Gaussianity that seems to be already ruled 
out by CMB constraints ( Komatsu et al. 2011), unless a strongly 
scale-dependent non-Gaussianity is invoked (as proposed by e.g. 
ILo Verde et al.ll2008l; ICavon et al.jl2010l ; iHovle et al.ll201 ll) . Alter- 
natively, ISaldi & Pettorind ( 1201 ih have shown that interacting DE 
models, which are generically characterized by a faster growth of 
density perturbations with respect to ACDM due to the presence of 
a fifth-force mediated by the DE scalar field, also predict a larger 
number density of massive halos at all redshifts that would result in 
a higher detection probability. 

Both these approaches, however, are faced by the additional 
problem of overpredicting the abundance of massive clusters at 
low redshift. In fact, while the detection of massive clusters at 
high redshifts seems to imply an anomalous growth of density per- 
turbations as compared to the predictions of the standard model, 
the observed number counts at low redshifts are found to be in 
very good agreement with t h e predictions of ACDM (see e.g. 
Reiprich & Boehringep I2OO2I; I vikhlinin et ahl |2009| ; iMantz et al.l 
201(]| ; lRozo et alj201 0). In other words, a non-Gaussian initial den- 



sity field or a faster growth of density fluctuations can determine 
an increased number density of massive halos at high redshifts, 
thereby alleviating the possible tension with the model arising from 
the unexpected detection of an exceeding number of such objects, 
but will then also necessarily imply the further growth of such mas- 
sive halos from these high redshifts to the present time, resulting in 
a clear discrepancy with present bounds on the cluster mass func- 
tion at low redshifts. 

In the present study we will show how this problem can 
be naturally solved by an interacting DE model with a suitable 
self-interaction potential different from the standard exponential 
or power-law potentials assumed in previous investigations of 
coupled DE (cDE) cosmologies. In particular, we will show 
how the existence of a global minimum in the potential, and the 
consequent "bounce" of the DE field on the cosmological constant 



barrier, is the key feature to address the existence of exceedingly 
massive clusters at high redshift without necessarily affecting the 
halo mass function at the present epoch. 

Although our proposed scenario does not benefit in general 
of the same scaling properties typical of the exponential and 
power-law potentials, it naturally provides a mechanism to en- 
hance the expected number of massive halos up to some high 
redshift and to subsequently reduce it again to the standard ACDM 
prediction at the present time, thereby accounting at the same 
time for the existence of massive clusters at early epochs and for 
the observed halo abundance at low redshift. We therefore regard 
the specific model presented in this work as a toy example of 
how the dynamical nature of a DE scalar field interacting with 
cold dark matter (CDM) particles can account for the presence 
of a high redshift peak in the deviation of the expected cluster 
number counts from the ACDM prediction. This feature could not 
arise in standard gravity models or in non-Gaussian cosmological 
scenarios, and would therefore represent a "smoking gun" for the 
dynamical nature of DE. In the present work we will illustrate 
in detail this peculiar behavior, also by means of large N-body 
simulations of structure formation. 

The paper is organized as follows. In Section |2] we briefly 
review the main features of standard cDE models for the back- 
ground evolution of the universe, and we highlight how bouncing 
cDE models deviate from this standard and well-known behavior. 
In Section[3]we investigate the evolution of linear density perturba- 
tions in the context of both standard and bouncing cDE scenarios, 
and we show how the latter can provide a possible explanation for 
a larger amplitude of density perturbations at high redshift without 
affecting the normalization of the linear matter power spectrum at 
the present time. This effect is further investigated for the nonlin- 
ear regime of structure formation in Section [4] where we show by 
means of large N-body simulations how the real halo mass function 
evolves in time in each of the models considered, and how bouncing 
cDE models can predict an enhanced number of massive clusters at 
high redshift without changing the halo abundance at the present 
epoch. Finally, in Section|5]we draw our conclusions. 



2 INTERACTING DARK ENERGY MODELS WITH 
EXPONENTIAL AND SUGRA POTENTIALS 

Interacting DE models involving a direct exchange of energy- 
momentum between a DE scalar field cj) and other matter com- 
ponents of the universe have been widely investigated in the 
past, in particula r for the case of a coupling with the CDM fluid 



past, in particula r tor tne case ot a couplmg witn tne LJJM nuia 
(Wetterich 199 ^ lAmendolall2000l 120041 ; IPettorino & Baccigalupi 
2008; Bal dii 120110) or w ith massive neutrinos JAmendola et al 



2008; Bal di et al.l uOllbl) . The effects of these interactions on 



observable quantities such a s the CMB an gular power spectrum 
JBean et al.i 2008; Xia 2wA iLa Vacca et al. 2009) , the Lyman-a 
forest JBaldi & Viell201(]n . the large scale structure of the universe 
I IBaldi & Pe ttoring20 1l|), or the s tructural properties of highly non- 
Unear objec ts (Baldi e7al.ll2010l ; ILi & Barrowiboi ll ; iBaldilbol id; 
iBaldi et al.l20Ilai') have been extensively studied by means of both 
analytic and numerical techniques. 

Most of the previous works on interacting DE models as- 
sume a monotonic function for the self-interaction potential of the 
DE scalar field, a s e.g. an exponen tial (Wetterich 1988) or an in- 
verse power law JRatra & Peebleslll988h potential, which deter- 
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Figure 1. Left: The trajectoiy of the DE scalar field </> along its self-interaction SUGRA potential during the whole expansion history of the universe. The field 
is set at rest in the minimum of the potential at 2 — > oo (blue square) and subsequently moves away from the minimum due to the effective force exerted by 
the coupling term in the Klein-Gordon equation (T). The field reaches the position marked by the green asterisk at ^CMB and keeps moving up the potential 
until the inversion point at zjnv ~ 6.8 (pink triangle) where it inverts its motion and rolls back down towards the minimum, reaching at z = the position 
marked by the red cross. Right: The same dynamical evolution described in the left panel: the upper plot shows the evolution in redshift of the scalar field <f> 
while the lower plot displays the corresponding evolution of the specific scalar velocity 4'/ H- 



mine the existence of attractor solutions that are almost indepen- 
dent on the initial conditions of the scalar field. These potentials 
do not possess local minima, and therefore define a unique direc- 
tion for the time evolution of the scalar field ij>. However, other 
possible functions, characterized in some cases by the presence of 
local or global minima, have been considered in the literature for 
the cas e of non-interac ting DE scenarios, as e.g. the SUGRA po- 
tential iBrax & Martial 1999,) arising naturally in supersymmetric 
theories of gravity. The existence of a local or a global minimum in 
the scalar field self-interaction potential allows the field to oscillate 
and invert its motion duri ng the expans i on his tory of the Universe, 
as recently suggested by IXarrant et al.l ( 120 111) . Such feature could 
have very significant consequences if a coupling to CDM is also 
present, as we will show in the present work. 

We will study the case of an interacting DE model with 
a SUGRA self-interaction potential and compare it with the 
standard ACDM cosmology and with coupled DE models with 
an exponential potential, with a particular focus on the different 
effects that these two types of potential induce on the evolution 
of linear and nonlinear density perturbations in the universe and 
on the predicted number density of massive clusters at different 
cosmic epochs. 

We consider a flat cosmological model described by the fol- 
lowing set of dynamic equations: 



-'«*+S ^ #-ife' 



pc + iHpc 

pb + 2>Hpb 

pr + '^Hpr 

3H^ 



2 pa(j} 



(1) 

(2) 

0, (3) 

0, (4) 

jp- (P'' + Pc + Pb + P4>) , (5) 



where the role of DE is played by the classical scalar field (f) and 
(c) 201 1 RAS, MNRAS 000,[T]tTT] 



where the subscripts 6, c, and r represent the baryonic, CDM and 
relativistic fluids, respectively. In Eqs. dllSt an overdot represents 
a derivative with respect to the cosmic time t, H = a/a, and 
Mpi = I/vSttG is the reduced Planck massQ In this work we 
will restrict our analysis to the case of a constant coupling /?c, al- 
though a direct dependence of the couplin g on the scalar field could 
also be considered (see e.g. iBaldiluOlld) . 

As already anticipated, besides the standard ACDM model 
we will consider two possible choices for the scalar field 
self- interaction p otential V{4>), namely an exponential poten- 
tial I Lucchin & M atarrese 1985; Wetterich 1988; Ferreira & Joyed 
ll998l) andaSUGRA potential JBrax & MartirJI 19991) , defined as: 



EXP : 1^(0) = ^e""^ , 
SUGRA: 1/(0) = yl(^!)""e"*'''^ 



(6) 

(7) 



where the scalar field is normalized in units of the reduced Planck 
mass. For the exponential models we will assume a positive cou- 
pling with values /?c = 0.1 and /3c = 0.15, while for the SUGRA 
model a negative coupling /3c = —0.15. Our proposed SUGRA 
cDE scenario therefore requires the same number of parameters as 
standard cDE models. 

We compute the background evolution of these models by 
numerically integrating the system of Eqs. ( IU5b with a trial-and- 
error procedure until the the desired values of the background 
cosmological parameters at z — are obtained. These correspond 
to the "WMAP7-only M aximum Likelihood" parameters derived 
bv lKomatsu et al.l ( I2OI li) . and are listed in Table [T] All the models 
considered in the present work with the relative parameters and 
normalizations are instead summarized in Table|2| 



^ Notice that in Eqs. JU5) and throughout the rest of this paper we assume 
the same convention on the definition of the coupling adopted bv lAmendola 
(2000) and Baldi (201 Ic). which differs by a factor y^3j2 from the defini- 
tion adopted in some of the literature. 
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Parameter 


Value 


Ho 


70.3kms-^Mpc-^ 


f^CDM 


0.226 


f^DE 


0.729 


0-8 


0.809 


Hi 


0.0451 


ris 


0.966 



Table 1. The cosmological parameters at 2 = assumed for all the models 
considered in the present study, consistent with the "WMAP7 only Maxi- 
mum Likelihood" results of .Komatsu et al.. (.201 1.) . 
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1.1 - 



1.0 r^ 



Since the main features of the cDE models with an exponential 
potential have been extensively discussed in the literature (see e.g. 
Wetterich 1995; Amendola 2000, 2004; Pettorino & Baccigalupil 
2008l : lBaldj201 id lai) we refer the interested reader to these publica- 
tions and we will mainly focus our attention here on the dynamics 
of the new proposed SUGRA cDE scenario which has never been 
presented before. The key feature of the SUGRA potential for the 
analysis carried out in this work is the presence of a global mini- 
mum at (j>m = \/Q' where the derivative of the potential in Eq. lUJ 
vanishes. For our model, we assume the scalar field to be at rest 
in the global minimum at very early times. This is a natural choice 
for the initial conditions since any previous dynamical evolution of 
the scalar field could be efficiently damped by Hubble friction in 
the radiation dominated epoch. 

In the absence of any coupling to CDM the global minimum 
of the SUGRA potential would be a stable critical point of the sys- 
tem, and the scalar field would therefore remain at rest at its initial 
location ^m during the whole expansion history of the universe, ef- 
fectively behaving as a cosmological constant. On the contrary, if a 
coupling to CDM is present, the Klein-Gordon equation (TJ for the 
scalar field at rest in the minimum reads: 



= 



^vl^^ifc 



(8) 



showing that a non-vanishing acceleration term due to the coupling 
acts on the scalar field as soon as the CDM density pc starts to 
play a significant role in the dynamical evolution of the universe. 
In particular, for our choice of the parameters in the SUGRA cDE 
model the negative coupling j3c = —0.15 will push the field away 
from its initial location (j)m in the direction of smaller field values, 
i.e. towards the power-law side of the SUGRA potential. As soon 
as the scalar field (j) leaves the local minimum with a negative ve- 
locity {(f) < 0) the other terms appearing in Eq. (TJ are restored 
and start to slow down the field, which will eventually stop and in- 
vert its motion moving again towards the minimum of the SUGRA 
potential. 

The details of this peculiar dynamic evolution depend cru- 
cially on the balance between the different parameters of the model, 
and in particular on the ratio between the slope of the SUGRA po- 
tential a and the coupling strength /3c. It is therefore necessary to 
tune these parameters quite accurately in order to obtain a viable so- 
lution for the background evolution of the universe, and the model 
therefore presents a comparable level of fine tuning as for the stan- 
dard ACDM cosmology. Nevertheless, as we will describe below, 
this scenario presents several appealing features and provides a self 
consistent way to address the issue of anomalous massive clusters 
at high redshift. It should therefore be considered as a toy example 
of the possible characteristic signatures of a dynamic DE field on 
structure formation. 

In the two panels of Fig. [T] we show the dynamic evolution 
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Figure 2. The time evolution of the CDM particle mass in the different 
cDE models. The different curves show the CDM particle mass in units 
of its present-day value. While the standard exponential cDE models show 
a monotonic behavior of the CDM mass, with a continuous flux of energy 
from CDM to DE. the SUGRA cDE model, due to the "bounce" of the scalar 
field and the inversion of its motion features a phase of CDM mass decrease 
followed by a subsequent increase. The trend is inverted in correspondence 
of the scalar field bounce, indicated in the figure by the vertical dotted line. 



of the scalar field (j) in the context of our proposed SUGRA cDE 
model. In the left panel we plot the evolution of the field along 
the SUGRA self-interaction potential at different times. The dot- 
ted curve shows the SUGRA potential for our choice of parameters 
while the solid line shows the trajectory of the field along the po- 
tential. The field is set to be at rest in the minimum of the potential 
<?!'m ~ \fci ~ 1.47 at early times (as shown by the blue square) and 
then evolves towards smaller field values reaching the position indi- 
cated by the green asterisk at the redshift of CMB, zcmb ~ 1100. 
The location of the scalar field at zcmb is particularly relevant for 
the analysis carried out in the present work, as will be extensively 
discussed below. After zcmb the field keeps then moving up the 
potential and eventually stops and inverts its motion at redshift Zinv 
(shown in Fig.[T]by the pink triangle) which for the specific model 
considered here is at Zinv ~ 6.8, and finally rolls back down to- 
wards the minimum, reaching at z = the position indicated by 
the red cross. It is important to notice here, for reasons that will be 
clarified later on, that the field values at zcmb and at z = are 
relatively close to each other 

The same dynamics are displayed in the right panel of Fig.[T] 
where the evolution of the scalar field (f> (upper plot) and of its spe- 
cific velocity (j>/H (lower plot) are shown as a function of redshift. 
Once more, the plot shows how the field starts at rest in the mini- 
mum of the SUGRA potential at high redshifts and then acquires a 
negative velocity due to the negative coupling term until the motion 
is eventually inverted and its velocity becomes positive at Zinv. The 
vertical dotted lines indicate the two redshifts zcmb and Zinv. 

The inversion of the motion of the scalar field during cos- 
mic evolution in our SUGRA cDE scenario, as opposed to stan- 
dard cDE models with monotonic potentials where the scalar field 
moves always in the same direction for the whole expansion his- 
tory of the universe, plays the crucial role in determining the effects 
on the growth of density perturbations that we will discuss in the 
next Section. In fact, as it is well kno wn (see e.g. Amendola 200C|; 
iBaldi et al j2O10l ; lTarrant et alj20I 11) . the dynamic evolution of the 
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Model 



Potential 



/3o /3i 



Scalar field 
normalization 



Potential 
normalization 



iy<^(2 = 0) As(zci./is) CTg(2 = 0) 



ACDM 


V{4>) = A 


- 


- 


- 


- 




A = 0.0219 


-1.0 


2.42 X 10-y 


0.809 


EXP0()2 


V{4') = Ae-°"'' 


0.08 


0.1 





<l,{z = 0) 


= 


A = 0.0218 


-0.995 


2.42 X 10-S 


0.875 


EXP0()3 


V{4>) = Ae-°"t' 


0.08 


0.15 





,l>{z = 0) 


= 


A = 0.0218 


-0.992 


2.42 X 10-S 


0.967 


SUGRA003 


V{</.) = A</.-"e'^''/2 


2.15 


-0.15 





4>{Z —>■ Cxd) : 


= VS 


A = 0.0202 


-0.901 


2.42 X 10-8 


0.806 



Table 2. The defining parameters and their normalization for the different cosmological models considered in the present work. All the models are consistent 
with the same amplitude of scalar perturbations As at ^CMB ~ 1100, but due to the different growth of perturbations the amplitude of the linear matter power 
spectrum at ^ = represented by erg is different from model to model. Notice that the SUGRA cDE cosmology is the only model to have the same amplitude 
of linear perturbations as the standard ACDM model both at zqmb and ai z = 0. 



DE scalar filed in the context of cDE models determines a transfer 
of energy-momentum between the DE and CDM fluids, according 
to the strength and sign of the coupling constant /3c. As one can 
clearly see from the CDM continuity equation (Eq. |2j, the trans- 
fer of energy between DE and CDM is directly proportional to the 
product of the coupling constant /?c and of the scalar field velocity 
(j). With our chosen convention, a positive combination /3c(^ > 
implies a flux of energy from the CDM to the DE component, while 
a negative combination J3c4> < G implies the opposite direction 
for the energy transfer. By integrating Eq. ^ and by assuming the 
conservation of the total CDM particle number one can therefore 
compute the time variation of the mass of CDM particles as: 



(^) 



rricfi ■ e 



(9) 



where mc,o denotes the value of the CDM particle mass at z = 0. 
It is therefore clear from Eq. l|9) that for the case of cDE models 
with a monotonic potential the CDM particle mass will be always 
decreasing or always increasing, while for the SUGRA cDE sce- 
nario the time derivative of the CDM particle mass ihc will change 
sign during cosmic evolution in correspondence to the inversion 
of motion of the DE scalar field, at Zinv A similar behavior has 
been recently describ ed also in the con text of Growing Neutrino 
quintessence models JBaldi et alJuOl Ibl) . This is shown in Fig. |2] 
where the evolution of the CDM particle mass for the ACDM, 
EXP002, EXP003, and SUGRA003 models (solid black, dotted 
green, dot-dashed blue, and dashed red lines, respectively) is plot- 
ted as a function of redshift. As one can see from the plot, while the 
mass is always decreasing in the standard EXP002 and EXP003 
cDE models, indicating a continuous flux of energy from the CDM 
to the DE, the mass in the SUGRA003 model first decreases (when 
< 0) and then increases again (when (j) > 0), such that the CDM 
particle mass at zcmb has almost the same value it has at z — 0, 
which is clearly not the case for a standard cDE scenario with a 
monotonic potential. This is obviously a consequence of the fact 
that for our specific choice of parameters and initial conditions the 
final location of the scalar field is relatively close, as mentioned 
above, to its location at the time of CMB. Such peculiar evolution 
of the CDM particle mass will have a significant impact on the 
growth of density perturbations, as we will discuss in the next 
Section. 

Before moving to consider the formation of cosmic structures 
in the context of our proposed SUGRA cDE model, we conclude 
the analysis of the background evolution by computing the Hub- 
ble function and the DE equation of state parameter, defined as 
w^ = {4P/2 - V{(t)))/{4)'^/2 + V{4))), for the four models con- 
sidered in our comparison. In Fig. |3] we show in the left panel 
the evolution of w^ with redshift for the EXP002, EXP003 and 
SUGRA003 models. As one can clearly see, the bounce of the 
scalar field at Zinv in the SUGRA cDE model, where the scalar 



field velocity cj) changes sign, is reflected in a bounce of the DE 
equation of state on the cosmological constant barrier at w^ = — 1. 
The following acceleration of the scalar field towards the minimum 
of the potential induces therefore a temporary increase of w^ up to 
w^ ^ — O.Satz ~ 2.5 until the Hubble friction starts decelerating 
again the scalar field motion, which ends up with an equation of 
state parameter w^ = —0.901 at 2 = 0. This non-trivial evolution 
of the equation of state parameter is not present for the standard 
exponential cDE models where w^ asymptotically approaches —1 
at low redshifts, and could therefore possibly strongly constrain the 
SUGRA cDE scenario. Its impact on the Hubble expansion, in fact, 
is quite significant at relatively low redshifts which could be probed 
by constraints on the expansion history coming from Supernovae 
la or from accurate m easurements of the Baryo n Acoustic Oscil- 
lations scale (see e.g. ISimpson & Bridlell2006l) from present and 
future wide surveys. This is shown in the right panel of Fig. [3] 
where the ratio of the Hubble function for the EXP002, EXP003 
and SUGRA003 cDE models (dotted green, dot-dashed blue and 
dashed red lines, respectively) to the standard ACDM case is plot- 
ted as a function of redshift. As one can see from the plot, the 
SUGRA cDE model determines a stronger difference from the stan- 
dard ACDM Hubble function with respect to an exponential cDE 
model with the same coupling strength at low redshifts, with a max- 
imum deviation of about 6% at z ~ 1, while at high redshifts the 
SUGRA model is found to be closer to ACDM as compared to the 
exponential case. 

This low-redshift feature of the SUGRA cDE might challenge 
our model by putting it in tension with observational constraints on 
the local expansion history. Nevertheless, a significant dynamical 
evolution of the DE scalar field at 2 < 2 is a necessary condition in 
order to account for the existence of anomalously massive objects 
at high redshifts and simultaneously recover the ACDM amplitude 
of density perturbations at 2 = 0, as we will discuss in full detail 
in Section[3] Furthermore, it is important to stress here that we are 
considering only one specific and somewhat extreme realization of 
the SUGRA cDE scenario, in order to enhance the effects under 
investigation, while a detailed parameter fitting would be required 
to select the set of parameters that best reproduce presently avail- 
able background observations. Ultimately, also the specific form of 
the self-interaction potential could be changed, since for our argu- 
ments the choice of the SUGRA potential is motivated mainly by 
its firm theoretical origin in the context of supersymmetric theo- 
ries of gravity, while also other functional forms of the potential 
could work equally well, provided they feature a global minimum. 
We are therefore proposing this specific model as a toy example of 
how a global or a local minimum in the DE self-interaction poten- 
tial could provide a mechanism to explain unexpected detections 
of massive clusters at high redshift, while a detailed survey of such 
models would be required in order to assess their viability for the 
background evolution of the universe. 
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Figure 3. Left: The evolution of the DE equation of state parameter ui^ as a function of redsliift. While the standard exponential cDE models show a monotonia 
evolution of w^ that asymptotically approaches the cosmological constant value of — 1 at low redshifts, the SUGRA cDE model features a "bounce" of the 
equation of state parameter on the cosmological constant barrier (sometimes also called the phantom divide) at Zi„v ~ 6.8, before decreasing again towards 
— 1 at lower redshifts, reaching at the present time a value ofw^(z = 0) = —0.901. Right: The ratio of the Hubble function H{z) over the standard ACDM 
case as a function of redshift. The non-trivial dynamics at low redshifts of the DE scalar field </< in the SUGRA cDE model imprints a specific signature on 
the Hubble function. For this model the deviation from ACDM is larger at very low redshifts as compared to standard exponential cDE scenarios, reaching a 
deviation of ~ 6% at 2: = 1, while at very high redshifts the deviation from ACDM is smaller with respect to a standard cDE model with identical coupling 
strength (EXP003). 



3 LINEAR PERTURBATIONS 

The evolution of linear density perturbations in the context of cDE 
cos mologies has been extensively studied by several authors (see 
e.g. Amendol j |2004| : IPettorino & Baccigalupil I2OO8I : iBaldi et al.l 
1201(1 ; lBaldil201 Ich to which we refer for the derivation of the main 
perturbations equations in the presence of a coupling between DE 
and CDM. In the present Section we recall the main modifications 
introduced by the coupling on the growth of density perturbations 
in the linear regime and we highlight how the dynamical evolution 
of a SUGRA cDE scenario extensively described above can sub- 
stantially alter this process and introduce noticeable features in the 
CDM growth factor. 

As a starting point, we consider the evolution equation for 
CDM and baryon density perturbations 5c,b = 5pc,b/pc,b in the 
Newtonian lim it and in th e prese nce of a DE-CDM interaction as 
derived by e.g. lAmendolal ( I2004I) . in Fourier space and in cosmic 
time: 



-2H 



Pc 



St = -2HSb + AtvG [pbSb + PcSc] 



5, + A-kG [pb5b + Pc5cTc] , (10) 

(11) 



The additional contribution appearing in the first term on the right 
hand side of Eq. JlOt is the extra friction a ssociated wit h momen- 
tum conservation in cDE models (see e.g. iBaldiluOllal for a dis- 
cussion on the effects of the friction term) and the factor Fc defined 



= 1+3/3.^ 



(12) 



includes the effect of the fifth-force mediated by the DE scalar field 
for CDM density perturbations. 

A few assumptions have to be made in order to derive 
Eqs. l llOlllt . besides the already mentioned Newtonian limit of 
small scales for the Fourier modes under consideration, i.e. A = 



aH/k <C 1. In particular, one has to assume the dimensionless 
mass of the scalar field cfi, defined as 



.2 I d'V 



H 



(13) 



to be small compared to the Fourier modes at the scales of interest, 
such that 



rh^A < 1. 



(14) 



If this condition is not fulfilled, the clustering of the DE scalar field 
4> might grow beyond the linear level at scales below 1/m^ and the 
fifth force of Eq. ( I12t would then acquire a Yukawa suppression 
factor given by: 



1 + 1 



Pl 



3 1 + m2A2 



(15) 



Since both the linear treatment of density perturbations discussed 
in the present Section and the non-linear N-body algorithm used for 
the analysis presented in the next Section are based on the assump- 
tion of a small scalar mass, and therefore on a term like J12b for 
the fifth-force implementation, it is important to clarify to which 
extent this assumption can be considered to hold. In Fig. |4] we 
plot the quantity m|A^ for several comoving wavenumbers k ^ 
1 , 0.1 , 0.01 , 0.001 /i/Mpc as a function of redshift. As one can 
see from the plot, the scalar mass can start playing a significant role 
at z < 10 only for scales close to the cosmic horizon, while for all 
scales below ~ 700 comoving Mpc//i (fc ~ 0.01 /i/Mpc) the in- 
fluence of a non-zero scalar mass is negligible and Eqs. il0lllll2l l 
can be safely used. For the N-body results discussed in Sec.|4]we 
will consider simulations with a box of 1 comoving Gpc//i aside, 
for which only the largest scales sampled by the initial power spec- 
trum might be marginally affected by our massless field approxi- 
mation, while we will concentrate on nonlinear structure formation 
processes occurring at much smaller scales. 
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Figure 4. The time evolution of tlie quantity mi A defining the validity 
of the linear approximation adopted both for the integration of the linear 
growth factor and for the N-body implementation used in the present work. 
As the Figure shows, fri^X^ gets close to unity only at extremely large 
scales (fc ~ 0.001 fc/Mpc). At the scales of interest in this work, the linear 
approximation assumed to derive Eqs. UQIIU is therefore fully justified. 
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Figure 5. The ratio of the linear growth factor D+ of the different cos- 
mologies considered in this work to the standard ACDM case. As the figure 
shows, all the models share the same perturbations amplitude at zcMB. 
while their subsequent growth significantly differs from model to model. 
The most noticeable feature of this evolution, which constitutes one of the 
central outcomes of this work, is that the SUGRA cDE model, differently 
from the standard exponential cDE scenarios, features a clear peak in the 
deviation of its growth factor from ACDM at ^inv ~ 6.8, while it recovers 
the same amplitude of linear fluctuations as ACDM at 2 = 0. 



Having clarified the range of validity of the assumptions made 
in deriving Eqs. dlOlllllJt we can now discuss one of the cen- 
tral results of the present work. By numerically solving the system 
of Eqs. l IllSllOllll at subhorizon scales {k ~ 0.1 fe/Mpc) we can 
compute the linear growth of density perturbations for all the dif- 
ferent cosmological models under investigation and compare it to 
the standard ACDM growth factor. As boundary conditions for our 
integration we impose that the ratio of baryonic to CDM perturba- 
tions at 2CMB takes the value 5b /5c ~ 3 .0 x 10~^ as com puted 
by running the Boltzmann code CAMB JLewis et alJuOOOf) for a 
ACDM cosmology and for the WMAP7 parameters adopted in the 
present study as listed in Table [T] The result of the integrations is 
shown in Fig. |5] where the ratio of the total growth factor D+ of 
each model over the standard ACDM case is plotted as a function 
of redshift. As one can see, all the models start from the same am- 
plitude of density perturbations at zcmb , but show a significantly 
different growth afterwards. 

The standard exponential cDE model EXP002 and EXP003 
(dotted green and dot-dashed blue lines, respectively) feature a 
faster growth during the whole expansion history of the universe, 
reaching at 2: = an amplitude of density perturbations exceeding 
the ACDM prediction by ~ 8% and ~ 20%, respectively. This 
is a very well known effect that has been widely discussed in 
the li t erature (see e.g. JAmendolal I200CI : IPettorino & Baccigalupil 
120081 : iBaldi et al.l l2010i) . and that arises as a combination of 
the modified background expansion of the universe and of the 
extra friction and fifth-force terms in Eq. JlOt that contribute to 
accelerate the growth of density perturbations. Such enhanced 
structure formation process has be en shown to account fo r an 
early formation of massive clusters JBaldi & PettorinduOl If) but 
the faster growth until the present time also necessarily implies 
a significantly higher normalization for the linear matter power 
spectrum at z = 0, with a value of erg = 0.875 for EXP002 and 
as = 0.976 for EXP003, as compared to the assumed a^ = 0.809 
for the ACDM cosmology. Such a high normalization at 2: = 
could be in contrast with present measurements of the matter 



power spectrum at large scales (see e.g. iMontesano et alJuOl ih 
and represents the main unaddressed problem when trying to 
explain the existence of high-z massive clusters in the context of 
standard cDE models. 

The central result of the present work is that such problem 
might be avoided by bouncing cDE models, and in particular by 
the SUGRA cDE scenario discussed in this paper. In Fig. [5] the 
dashed red line shows the evolution of the ratio D+/D^ for 

our SUGRA003 cosmology. While at early times the growth pro- 
ceeds in a similar fashion as for the standard exponential cDE 
model EXP003, this trend is suddenly stopped and inverted when 
the scalar field (j> changes its direction of motion at Zinv After this 
redshift the growth of density fluctuations is slowed down below 
the rate of the ACDM cosmology such that the discrepancy be- 
tween the two models is progressively reduced and the amplitude 
of perturbations reaches again the ACDM value at 2 = 0. This 
implies that the linear power spectrum normalization of the two 
models will be the same both at 2cmb and at 2 ~ 0, while a sig- 
nificantly larger amplitude is reached in the SUGRA model at an 
intermediate redshift 2inv 

This peculiar evolution of the matter growth factor could not 
be realized in the context of standard gravity models or by a mod- 
ification of the statistical properties of the initial conditions for the 
density fluctuations as for the case of non-Gaussian cosmological 
scenarios. It is therefore a quite specific footprint of a non trivial 
dynamical behavior of the DE scalar field at relatively late cosmo- 
logical epochs. As we will show in the next Section, this feature of 
the SUGRA cDE model can account for a higher number density 
of massive clusters at early times while not significantly affecting 
the halo abundance at 2 = 0. 

Furthermore, since the different growth of structures induced 
by the DE interaction is expected to affect in a non- trivial way also 
the peculiar velocity field (see e.g. iLee & Baldi 2011), a detailed 
investigation of the impact of cDE models on redshift-space dis- 
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tortions in the large-scale distribution of galaxies mi giit provid e an 
independent way to test the cDE scenario cMarulli et alj201 ih . 



4 NON-LINEAR STRUCTURE FORMATION 

In order to test directly the effects that the SUGRA cDE scenario 
has on the expected number density of massive halos at different 
cosmic epochs we need to rely on the results of large N-body simu- 
lations which should include in their numerical treatment a self- 
consistent implementation of the cDE effects described in Sec- 
tion [2] To this end we make use of the pu blicly availab le halo 
catalogs of the CoDECS simulations suitcl ( lBaldill20libh which 
includes the very same EXP002, EXP003 and SUGRA003 mod- 
els investigated in the present work. The CoDECS simulations have 
been performed using the modified version bv lBaldi e t al. ( 201(|) of 
the parallel Tree-PM N-body code GADGET-2 ( ISpringel.2005.) , that 
implements all the characteristic features of cDE models described 
above. We refer to the CoDECS website for an extensive description 
of the numerical setup of the simulations, and we recall here only 
their main features. 

The ACDM-L, EXP002-L, EXP003-L, and SUGRA003-L 
simulations of the CoDECS database that are considered in the 
present work consist of a cosmological box with a size of 1 co- 
moving Gpc//i and periodic boundary conditions, filled with 1024^ 
CDM particles with amass of mc = 5.84 x 10^" Mq/H, and with 
1024^ baryonic particles with a mass of mt = 1.17x 10^" Mq//i. 

Initial conditio ns are gene rated by perturbing a homogeneous 
glass distribution ( Whitdl 19941) in order to set up a random-phase 
realization of the initial linear matter power spectrum which is 
taken to be the spectrum c omputed by the publicly available code 
CAMB j Lewis et aLll2000l , www.camb.info) at the starting redshift 
of the simulations Zi = 99 for the WMAP7 parameters specified 
in Table [T] A common normalization of the linear perturbations 
amplitude at zcmb consistent with present WMAP7 constraints is 
imposed to all the simulations by properly rescaling the particles 
displacement between zcmb and Zi with the specific growth factors 
numerically computed for each cosmology. This procedure makes 
the Codecs runs more suitable to realistically quantify the impact 
of cDE o n the halo mass function w ith respect to previous attempts 
(as e.g. bv lBaldi & Pettorincl201in where a common normalization 
at the starting redshift of the simulations was assumed. A more ac- 
curate quantitative determination of the number counts deviation 
from ACDM can therefore be expected by using the CoDECS halo 
catalogs. These are obtained by running a Friends-of-Friends (FoF) 
algorithm with linking length ^ = 0.2 x J, where d is the mean 
interparticle separation. Since the CoDECS runs include both CDM 
and baryonic particles, the halo catalogs are compiled by running 
the FoF algorithm on the CDM particles as primary tracers and 
then linking the baryonic particles to the FoF group of their closest 
CDM neighbor. It is also important to stress here that all the simu- 
lations are started with the same initial linear transfer function and 
with the same random seed for the realization of the matter power 
spectrum in the initial conditions. 

Having at hand all the CoDECS FoF halo catalogs, it is 
possible to compute the cumulative halo mass function for each of 
the models under discussion at different redshifts and to compare 
their relative evolution. This is shown in the three panels of Fig.|6] 



^ The CoDECS simulations are publicly available at the URL: 
www.marcobaldi.it/research/CoDECS 
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Figure 6. The halo mass function as extracted from the halo catalogs of 
the CoDECS database, at 2 = 2.5 (top), 2 = 1.6 (middle), and 2 = 
(bottom). For each panel, the upper plot shows the cumulative halo mass 
function in the mass range covered by the sample at that redshift, while 
the lower plot shows the enhancement of the expected number counts with 
respect to ACDM. While all the cDE models predict a larger number of 
massive halos at 2 > 0, our proposed SUGRA cDE scenario is the only one 
that recovers the ACDM expectations at the present epoch, thereby avoiding 
tension with available bounds on the cluster abundance at 2 = 0. 
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for three choices of redshift, namely z — 2.5, z — 1.6 and z = 0. 
For convenience, the mass functions are plotted by binning halos 
into 15 logarithmically equispaced mass bins over the whole mass 
range covered by the sample at each redshift. 

As one can clearly see from the Figure, at redshifts as high 
as z — 2.5 all the cDE models show a larger number density 
of halos with respect to ACDM over the whole mass range cov- 
ered by the catalogs. Th is effect was already shown in the past by 
iBaldi & Pettorind BoTlh but the results of the CoDECS simulations 
significantly extend its validity to higher masses, besides provid- 
ing a more realistic quantitative determination of the overall effect 
thanks to the common CMB normalization of the linear amplitude 
in all the models. The increase in the halo number density, dis- 
played in the bottom plot of each panel as the ratio of the halo mass 
function of each model over the ACDM case, shows a clear mass 
dependence, being more significant at large masses and suggesting 
a possible further extrapolation of these results to more massive 
halos than the ones found in the CoDECS catalogs. In particular, 
the EXP003 model shows an increase of almost a factor of 20 at 
M ~ 1.0 X 10^* Mq/Zi while at this redshift the EXP002 and 
SUGRA003 models determine a comparable increase of a factor 
~ 4 — 5 at the same mass. An accurate fit of the mass fun ction, 
that we defer to an upcoming publication JCui et afllin prepQ . will 
clearly be required in order to determine a reliable extrapolation of 
these numerical results to larger masses and to investigate in detail 
to which extent the deviation in the number counts is degenerate 
with the evolution of the linear growth factor. Nevertheless, a rough 
estimate suggests that both the EXP002 and the SUGRA003 cos- 
mologies might reach at z — 2.5 an enhancement factor of ~ 10 
with respect to the ACDM halo mass function at A/ ;^ 3 — 6 x 10^'' 
Me/h. 

At later times, as one can see from the central panel referring 
to z — 1.6, the situation is already significantly different. While the 
EXP003 model still provides the strongest effect, reaching a factor 
of ^ 10 in the increase of the halo mass function at A/ ~ 2 x 10^^ 
Mo//i, the EXP002 and SUGRA003 models now show a signifi- 
cantly different amplitude of the number counts enhancement: the 
former reaches an increase of a factor ~ 6 at Af ~ 2.5 x 10^^ 
M0//i, while a factor of ~ 3 is attained by the latter at the same 
halo mass. Also in this case, the enhancement shows a clear 
mass trend in all the models, and a larger effect can be expected 
for halo masses beyond the range covered by the CoDECS catalogs. 

The central result of the present work is shown by the direct 
comparison of the top and central panels of Fig. [6] just discussed 
with the bottom panel, that displays the halo mass functions of the 
different models at 2 = 0. While both the standard exponential 
cDE models EXP002 and EXP003 still show at the present time a 
similar enhancement of the halo number density as found at higher 
redshifts, with an increase at A^ « 3.0 x IO^^Mq/Zi of a factor of 
~ 4 and ~ 15, respectively, the bouncing cDE model SUGRA003 
does not show any excess in the halo number density over the whole 
mass range of the sample, and is therefore found to be perfectly 
consistent with the ACDM predictions on the cluster abundance at 
z — 0. This is a very remarkable result, since it shows how the 
dynamic degree of freedom associated with a cDE scalar field 
can provide, for suitable self-interaction potentials, a natural mech- 
anism to explain an excess of massive clusters at high redshift as 
compared to the ACDM expectations without affecting the very 
tightly constrained cluster mass function at the present time. Such 
behavior could not easily arise in the context of standard gravity 



theories or non-Gaussian cosmological scenarios, and would there- 
fore represent a clear indication of a dynamical nature of DE. 

Further detections of anomalously massive clusters at high 
redshift could therefore be considered as a "smoking gun" for the 
existence of a dynamic degree of freedom in the dark sector. 



5 CONCLUSIONS 

In the present paper we have proposed a new model of interacting 
dark energy based on a SUGRA self-interaction potential and on 
a constant coupling between dark energy and CDM particles. The 
most relevant feature of the SUGRA self-interaction potential for 
the analysis carried out in this work, is the presence of a global min- 
imum that allows the dark energy scalar field to oscillate, thereby 
changing its direction of motion during the expansion history of 
the universe. Such feature makes our SUGRA coupled dark energy 
scenario significantly different from previously proposed models 
of dark energy interactions, as it allows an inversion of the energy 
flow between the two interacting components at relatively recent 
cosmological epochs. 

We have therefore compared the cosmological evolution of 
standard coupled dark energy models based on an exponential 
self-interaction potential to our proposed SUGRA cDE scenario, 
both concerning the background dynamics and the evolution of 
linear and nonlinear perturbations. For the latter task, we have also 
made use of large N-body simulations results through the publicly 
available CoDECS halo catalogs that include the specific models of 
dark energy interaction under investigation. 

Already at the background level, the SUGRA cDE scenario 
shows a series of interesting features. First of all, the inversion of 
the direction of motion of the dark energy scalar field happens in 
our specific model at Zinv ~ 6.8 which is a redshift already rele- 
vant for astrophysical processes and for the early stages of nonlin- 
ear structure formation. This inversion of motion corresponds to a 
"bounce" of the dark energy equation of state parameter w^ on the 
cosmological constant barrier w^ — —1, and imprints a specific 
pattern in the low-redshift evolution of the Hubble function that 
might provide a way to directly constrain the model. Furthermore, 
the inversion of the energy flow between dark energy and CDM 
particle during cosmic evolution implies, for the specific model pre- 
sented here, that the CDM particles mass has very similar values at 
the redshift of the CMB and at the present time, thereby avoiding 
any mismatch of the cosmological parameters between zcmb and 
z = 0. 

The most significant result of the present work, however, 
concerns the evolution of perturbations and the formation of 
linear and nonlinear structures in the context of the SUGRA cDE 
scenario, as this is shown to provide appealing features to account 
for possible anomalous detections of massive clusters at high 
redshifts. 

At the linear level, we have demonstrated for the first time 
by numerically solving the linear perturbations equations that 
a bouncing coupled dark energy model, and in particular the 
SUGRA cDE scenario presented here, can determine the same 
amplitude of linear perturbations as a corresponding ACDM 
cosmology both at CMB and at the present time, while featuring 
a significant deviation from ACDM at intermediate redshifts. In 
particular, both the amplitude of scalar perturbations at CMB and 
the value of cts at 2; = are the same between the standard ACDM 
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model and our proposed SUGRA cDE scenario, while the latter 
features a significantly larger amplitude of density perturbations at 
intermediate redshifts, with a peak in correspondence of the dark 
energy "bounce" at Zinv Therefore, the SUGRA cDE model can 
be simultaneously consistent with CMB and local constraints on 
the perturbations amplitude, while allowing for significant freedom 
at intermediate epochs. 

We have then studied the evolution of the halo mass function 
as predicted by large N-body simulations that include all the cha- 
racteristic features of the different interacting dark energy models 
considered in this work. Consistently with what found in previous 
studies, at high redshifts (z ~ 2.5) all the coupled dark energy 
models show a larger number of halos with respect to ACDM over 
the whole mass range covered by our numerical sample. The en- 
hancement in the halo number density reaches a factor of ~ 20 
at M ~ 1.0 X 10^'* Mq/Zj for the most extreme standard cou- 
pled dark energy scenario, while the SUGRA cDE model shows 
an enhancement of a factor ~ 4 — 5 at the same mass. Further- 
more, the enhancement has a clear mass dependence and increases 
towards larger masses, suggesting that the effect might be signif- 
icantly larger at masses not covered by the halo sample of our 
simulations. The subsequent evolution of the halo mass function 
clearly displays the fundamental difference between the standard 
coupled dark energy models based on a monotonic self-interaction 
potential and the bouncing coupled dark energy scenarios, as the 
SUGRA cDE proposed here. At 2; ~ 1.6, in fact, all the models 
still feature a significant excess of massive halos as compared to 
ACDM, although the SUGRA cDE model has reduced its enhance- 
ment factor much more significantly than the other coupled dark 
energy models. Finally, at 2 = the SUGRA cDE model fully 
recovers the standard ACDM mass function over the whole mass 
range of our sample, while the other coupled dark energy models 
still show a large excess of massive halos, especially in the range of 
very massive galaxy clusters, being therefore potentially in tension 
with available constraint on the cluster number counts. 

This result, which represents the main outcome of the present 
paper, demonstrates for the first time by means of large and 
fully self-consistent N-body simulations that coupled dark energy 
models with a suitable choice of the self-interaction potential can 
simultaneously account for the detection of anomalously massive 
clusters at high redshifts and for the observed halo mass function 
at low redshifts. Such peculiar evolution could not arise in other 
types of models that have been recently invoked as a possible 
explanation of the unexpected detection of massive clusters at high 
redshifts, as e.g. non-Gaussian cosmological scenarios, where an 
enhanced number density of halos at high redshifts necessarily 
implies a corresponding enhancement at low redshifts. 

To conclude, we have studied a new class of interacting dark 
energy cosmologies characterized by a "bounce" of the dark en- 
ergy scalar field that is allowed to invert its direction of motion 
during the cosmic expansion. We have shown by means of linear 
and nonlinear numerical treatments that this new class of models 
could be simultaneously consistent with observational constraints 
at CMB and at the present epoch, while allowing for significant 
deviations from the standard ACDM scenario at intermediate red- 
shifts. In particular, we have shown this class of models to possess 
the (so far) unique feature of simultaneously accounting for unex- 
pected detections of very massive clusters at high redshifts and for 
the standard cluster abundance at the present time. Such behavior 
is a direct consequence of a non-trivial dynamics of dark energy 



at relatively recent cosmological epochs, and therefore represents 
a specific observational signature of a possible dynamical origin of 
the accelerated expansion of the Universe. 
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